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ABSTRACT 


Householder's  method  reduces  a  general  matrix  to  an  upper  Hessenberg 
form  using  elementary  Hermitian  orthogonal  matrices.   For  a  matrix  of  order  n, 
there  are  n-2  steps  in  this  reduction.   An  ASK  subroutine  for  implementing 
this  method  on  ILLIAC  IV,  for  n  <  bh,    is  described  in  this  document.   A 
storage  scheme  for  maximizing  PE  efficiency  is  also  considered. 
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INTRODUCTION 


An  ASK  program   is  written  in  subroutine  form  for  reducing  a  real 
(n  x  n)  matrix,    n  <  Gh,    to  an  upper  Hessenberg  form  by  Householder's  method. 
This  method  is  one  of  the  most   efficient   and   stable  methods  to  reduce   a  gen- 
eral matrix  to   an  upper  Hessenberg  form.      The   ILLIAC   IV  computer  time  required 
for  reducing  a  real  6U  x  6k  matrix   into  this   form   is   approximately  200  milli- 
seconds . 

The  necessary  information  for  using  this   subroutine   is  given   in 
Section  2.      The  test   result  for   a  real  5x5  matrix   is  given  in  Section  3* 
In   Section  k  the  outline  of  the  Householder's  method   is  given.      A  new  modified 
skew  storage   scheme  which    is  used   in  this  program   is   described   in  Section  5- 
In  Section  6,    some   considerations   are  given  on   implementation  of  this   algorithm 
on  the   ILLIAC   IV. 


2 .   USAGE 

This  program  assumes  that  the  matrix  to  be  reduced  is  stored  in 
the  primary  memory  in  the  straight  storage  scheme.   This  means  that  each  row 
is  stored  across  PE's  starting  with  PEO.   The  reduced  matrix  replaces  the 
original  matrix  and  is  also  in  the  straight  storage  form. 

2.1  Calling  sequence 

Calling  sequence  for  this  subroutine  is: 
CALL   HSBG(N,A); 
where  A  designates  the  first  row  of  the  matrix  and  is  given  as 
A:  DATA     V  aQ1,  ...  ,  a  ^  (Q.O)M, 


a10'  all'  '•*  ' 


al,N-l'  (°*0^ 


Vi,c  Vi,i'  •'••  '  Vi,m'  (0-0)M* 


where     M  =  6^  -  N. 

N  is  the  size  of  the  matrix  declared  as 

N:   EQU    6k   j 

or        DEFINE  N  =  6k   #  #  ; 
or  given  as  an  integer,  i.e.,  6k. 

2.2  Core  storage  used  by  this  subroutine 

Besides  the  space  used  for  storing  instructions  which  occupies  226 
words,  this  subroutine  uses  2  rows  of  PE  memory  for  storing  PE  numbers  j's 
and  modified  PE  numbers  j's  (see  Section  5).   PE  numbers  are  stored  in  the 
row  named  PEN  and  j's  are  stored  in  the  row  PENBR.   For  instance  PENBR  con- 
tains 0,  32,  16,  k8,    .  .  .  at  PE  0,  1,  2,  3,  ...  ,  since  0  =  0,  1  =  32, 
2  =  16,  3  =  ^-8  and  so  forth.  Four  rows  of  PE  memory  are  used  for  storing 
intermediate  results.  ADBO  ~  32  are  also  used. 


2. 3  Subroutine  called  within  this  program 

The  square  root  subroutine  is  called  in  this  program.   This  sub- 
routine is  declared  as  an  external  routine,  EXTERNAL  SQjRT.   It  is  assumed 
that  this  subroutine  is  available  to  the  users  as  a  separately  compiled  pro- 
gram, and  computes  the  square  roots  of  the  floating  point  numbers  placed  in 
registers  A  whose  mode  bits  are  set.   This  following  linkage  is  used 
for  calling  the  subroutine: 

CLC  (3)  ; 

SLIT  (3)  SQRT  ; 

EXCHL  (3)  $1CR  ; 

The  resultant  square  root  numbers  are  assumed  to  be  placed  in 
registers  A  at  the  return  from  the  subroutine. 


3-   EXAMPLE 


As  a  test  of  this  program,  the  matrix 


A  = 


0.32 

0.27 

O.23 

0.32 

o.k 

0.25 

0.03 

O.71 

0.21 

0.17 

0.1+3 

0-73 

0.13 

0.37 

0.85 

0.53 

0.25 

0.51 

0.62 

0.16 

0.16 

O.65 

0.1+6 

0.56 

0.32 

was   selected.      Table   1  gives   the  result  obtained  on  the  ASK  simulator  on  the 
B5500.      Because  of  the  very  slow   computation  speed  of  the   SSK  simulator   (ratio 
of  speed  to   ILLIAC   IV  is  1  to  10   ),   we  were   contented  here  with  a  matrix  of 
5x5  for  testing. 


Table  1 


0.32  -0.537^5772^10^     0.3089124l3802  -O.O57085503293   -0.025564833528 

■  O.7I+I+2I+1+583I+5I+  1. 1+597361+ li+5 15  -O.U96797706958  0.019l+79l+605i+6     0 .  189676107I+29 

0.0  -0.678586766101    0.117102691998  0.019392163333  -0.21+8301539^86 

0.0  0.0                         O.66152I+1+86138  -0. 3^361+093297  -O.2712508UI+712 

0.0  0.0                     0.0  -0.1+15155981387  -0.1321+75013216 


k.      HOUSEHOLDER'S  METHOD 

For  an  W  x  N  matrix,  the  Householder's  method  consists  of  N  -  2 
steps  of  reduction.   In  the  r  th  reduction,  where  r  =  0,  1,  2,  ...  N-3;  the 
zeros  are  introduced  in  the  r  +  lth  column  without  destroying  the  zeroes  in- 
troduced in  the  previous  steps.   The  configuration  of  the  matrix  immediately 
before  the  r  th  step  is  illustrated  below  for  the  case  N  =  ^ ,    r  =  2. 


r-1 


X  X  X  X  X  X  X 

X  X  X  X  X  X  X 

0  X  X  X  X  X  X 

0  0  x  x  x  x  x 

0  0  x  x  x  x  x 

0  0  x  x  x  x  x 

0  0  x  x  x  x  x 


Hr-1 

c..i 

o    : 

r-l         r-1 

The  matrix  H  ,   is  of  order  (r+1)  x  (r+l)  and  in  Hessenberg  form,  b 

r-1  r-1 

is  a  column  matrix. 

The  matrix  P   which  transforms  A  _   to  A   is  introduced  as 
r  r-1       r 


0 
0, 


I 

.      0 

0 

I  - 

T 
2  w  w 
r  r 

where  I  is  an  (r+l)  x  (r+l)  identity  matrix. 
Then 


A  =  P  A   P 
r    r  r-1  r 


H    1 

r-1 

r-1  r 

o  :c 

.    r 

Q  B     ,Q 
^r  r-1  °r 

X 

X 

X 

.     X 

X 

X 

X 

X 

X 

X 

.     X 

X 

X 

X 

0 

X 

X 

.     X 

X 

X 

X 

0 

0 

X 

.     X 

X 

X 

X 

0 

0 

0 

.     X 

X 

X 

X 

0 

0 

0 

.     X 

X 

X 

X 

0 

0 

0 

.     X 

X 

X 

X 

where  C   (=  Q  b   )  is  null  except  for  its  first  element.   Therefore,  we 

have 

P=I-2wwT=I-u  uT/2K2    such  that 
r         r  r        r  r'   r 


u.    =0   for   i  =0,  ....  r 

i,r  '    ' 


Ur+l,r  "  ar+l,r  +    S^  (l) 

u.  =  a.     for  i  =  r+2,  ....  N-l 
i,r      i,r 

ajad  W-l    0  \ 

K  (  E  a,  J  (2) 


r 


-K2 


r+1   ' 

=  S2  T  a  .   S 
r    r+l,r  r 


In  the  equation  defining  u       the  sign  associated  with  S  is  taken  to  be 

r+±,  r  T 

that   of  a  so  that   cancellation  cannot  take  place   in  the  addition 

r+l,r 

a     .        +  S      .      The  new  value  of     a     .  is  T  S    . 

r+1,  r         r  r+1,  r  r 

In  this  program,    the  r  th   step  which   computes  A     from  A     ,    is 

carried  out   in  the   following  manner. 

2 

I.  S     and  2K       are  computed  by    (2), 

II.  the  vector  u     as  defined  by    (l)    is   formed, 

III.  the  vector  v     =  u   /2Iv  is   computed, 

T     T 

IV.  p   is  computed  by  the  relation  p  =  u  A  ,, 

V.  q  =  A  nu   is  computed. 
-±r    r_2_  r 

VI.  the  matrix  A   is  computed  as 

m  m 

A   =  A  -,-vp  -  (q  -  OL    u  )u 
r    r-1    r-^r   x^r    r  r  r 

where 

a     =  p  v  . 
r   rr  r 


5.   MODIFIED  SKEW  STORAGE  SCHEME 

Given  any  positive  integer  j  such  that  0  <  j  <  63,  this  number  can 

S      k         '     3      2      1 
be  expressed  in  the  binary  form  as  j  =  2yj,-  +  2  j,  +  2  J3  +2^2+2^l+J'n 

for  jQ,  j^  j2,  j  ,  3^,  j'5  =  0  or  1. 

We  define  J  as  J  =  25jQ  +  2  ^  +  23J2  +  22^  +   21^  +  j  .   It  is 

obvious  that  0  <  j  <  63  for  any  0  <  j  <  63  and  j  =  j '  if  and  only  if  j=j'. 

In  modified  skew  storage,  the  elements  in  each  row  of  a  matrix  are 
first  reordered  so  that  the  jth  element  in  the  row  is  mapped  into  the  j  th 
position.   This  matrix  is  then  skewed  into  an  ordinary  skew  storage.   Fig.  1 
illustrates  this  storage  scheme  for  16  PE's  instead  of  the  actual  Gh   PE's. 
The  mapping  of  an  element  a. .   into  PE  memory  is  expressed  as 

a..  =  loc  1,    PE  (j  +  i)  (mod  6h) . 

This  storage  scheme  will  allow  simultaneous  access  to  each  element 
of  any  2   consecutive  row  vectors  or  column  vectors  whose  order  is  less  than 
64/2  for  k  =  I,  2.,    3>  4;  5-  For  example,  four  consecutive  rows  or  columns 
can  be  accessed  simultaneously  if  the  order  of  the  matrix  under  consideration 
is  less  than  16. 


8 


UA 

H 

c— 

d 

OO 

OO 

H 

UA 

ON 

H 

J" 

rH 

VO 

o 

H 

OJ 

?l 

-4" 

CO 

o 

H 

»\ 

•s 

•\ 

•A 

■\ 

•N 

•\ 

•\ 

■s 

•n 

•N 

•\ 

»\ 

•s 

■^ 

O 

H 

OJ 

OO 

-=f 

LTN 

VO 

t— 

CO 

ON 

o 

H 

r^ 

OO 

J" 

UA 

rH 

i-i 

H 

H 

H 

H 

OO 

-d- 

O 

OJ 

O 

UA 

E*- 

H 

OO 

H 

LT\ 

ON 

rH 

H 

VO 

H 

OJ 

H 

-d 

CO 

H 

H 

•N 

*\ 

•\ 

". 

•N 

•s 

»\ 

•A 

•N 

•s 

•s 

•s 

»\ 

•A 

•\ 

*\ 

O 

H 

OJ 

OO 

-cf 

LTN 

VO 

c- 

CO 

ON 

o 

H 

H 
H 

OJ 
H 

ro 
H 

J- 

rH 

UA 
H 

H 

OO 

J- 

o 

OJ 

LTN 

H 

ro 

rH 

LTN 

OA 

H 

H 

VO 

H 

OJ 

H 

-d- 

CO 

O 

H 

t- 

*\ 

•N 

»\ 

•\ 

•A 

»\ 

»\ 

•A 

•\ 

"    *\ 

•N 

•\ 

•\ 

>s 

*\ 

»\ 

O 

rH 

OJ 

OO 

J" 

LTN 

vO 

D~- 

CO 

ON 

O 
H 

rH 

H 

OJ 
rH 

OO 
rH 

-d 
H 

UA 
H 

ro 

_* 

o 

OJ 

LTN 

rH 

ro 

H 

LT\ 

o\ 

H 

rH 

VO 

H 

OJ 

H 

-rj- 

CO 

o 

H 

t- 

H 

*\ 

»s 

*s 

■N 

•\ 

O 

H 

OJ 

OO 

J- 

UA 

VO 

t— 

CO 

CA 

O 

rH 

H 

OJ 
H 

OO 

H 

-4- 
H 

UA 
H 

ro 

-d 

O 

rH1 

LPs 

rH 

rH 

LTN 

ON 

H 

H 

VD 

rH 

OJ 

-d 

CO 

O 

H 

t- 

H 

ro 

o 

H 

OJ 

OO 

-=f 

LTN 

VO 

t— 

CO 

ON 

o 

H 

H 

rH 

OJ 
H 

OO 

H 

-=Sr 

rH 

UA 

H 

-=r 

O 

^ 

LTN 

H 

ro 

LTN 

0\ 

rH 

H 

vO 

H 

OJ 

J- 

CO 

O 

rH 

t- 

rH 

ro 

rH 

*\ 

•N 

*\ 

•s 

•a 

*\ 

»\ 

•N 

•\ 

•N 

•N 

•N 

»\ 

•\ 

•S 

•N 

o 

H 

OJ 

OO 

_=r 

LTN 

VO 

C- 

CO 

ON 

o 

H 

H 
H 

OJ 
H 

OO 

H 

-d- 
H 

UA 

H 

J- 

o 

OJ 

LTN 

H 

ro 

ON 

H 

rH 

vO 

H 

OJ 

rH 

-=r 

CO 

o 

H 

f- 

H 

OO 

H 

UA 

»\ 

•A 

•s 

»\ 

•N 

•s 

•N 

•N 

•N 

•N 

»\ 

•\ 

*\ 

•A 

•A 

*\ 

o 

H 

OJ 

OO 

-=h 

LTN 

vO 

t- 

CO 

ON 

O 
H 

rH 

H 

OJ 

rH 

OO 
H 

-d 

H 

UA 

H 

J- 

o 

OJ 

LTN 

H 

OO 

rH 

H 

vO 

H 

OJ 

H 

-=J- 

co 

o 

H 

[— 

H 

OO 

rH 

LTN 

OA 

•A 

•\ 

•S 

•\ 

•A 

»\ 

•A 

•s 

»s 

■S 

•s 

»\ 

•N 

•s 

•\ 

•N 

o 

H 

OJ 

OO 

-J" 

LT\ 

VO 

t- 

CO 

ON 

o 
H 

rH 
rH 

OJ 
H 

OO 

H 

rH 

UA 
rH 

-=f 

o 

r^ 

LTN 

rH 

OO 

rH 

VO 

rH 

OJ 

J" 

CO 

o 

H 

t~- 

rH 

OO 

rH 

UA 

ON 

H 

•\ 

•A 

»\. 

*\ 

•\ 

*\ 

•A 

•N 

•s 

»s 

•\ 

•N 

•\ 

•\ 

»\ 

»\ 

o 

H 

OJ 

OO 

-=t" 

LTN 

VO 

t*- 

CO 

On 

o 

rH 

r^ 

OO 

-d- 

UA 

H 

H 

H 

rH 

H 

O 

rH1 

LTN 

rH 

OO 

J- 

VO 

rH 

OJ 

-3" 

CO 

O 

H 

t- 

H 

OO 

H 

LT\ 

ON 

H 

H 

O 

rH 

OJ 

OO 

J" 

LT\ 

VO 

t~- 

co 

OA 

O 

rH 

rH 
rH 

OJ 
H 

OO 

■H 

H 

UA 
H 

O 

r^ 

LTN 

H 

OO 

J- 

H 

OJ 

-J- 

co 

O 

H 

t- 

H 

OO 

H 

LTN 

ON 

rH 

rH 

VO 

*\ 

•N 

^\ 

*\ 

•N, 

•N 

*\ 

•N 

■N 

•N 

•N 

•N 

•s 

■\ 

•\ 

•\ 

O 

H 

OJ 

OO 

-=f 

LTN 

vO 

t- 

CO 

ON 

O 
H 

H 

H 

OJ 

rH 

OO 

H 

-d 

rH 

UA 

H 

OJ 

UA 

rH 

OO 

J- 

O 

OJ 

rH 

-3- 

CO 

o 

rH 

t~- 

rH 

OO 

H 

UA 

ON 

rH 

rH 

vO 

H 

»s 

»\ 

*\ 

■s 

•\ 

•\ 

*\ 

•N 

*s 

•s 

•\ 

•N 

•N 

•\ 

•N 

»\ 

o 

H 

OJ 

OO 

_=f 

LTN 

VO 

t- 

CO 

ON 

O 
H 

rH 
rH 

OJ 

rH 

OO 
rH 

rH 

UA 

H 

OJ 

LT\ 

H 

OO 

_* 

O 

H 

-d 

00 

o 

rH 

c— 

H 

OO 

rH 

LTN 

ON 

rH 

H 

VO 

H 

OJ 

■N 

•\ 

*\ 

•> 

•N 

•A 

»v 

•\ 

•\ 

*\ 

»\ 

*\ 

•\ 

»s 

»\ 

*\ 

O 

H 

OJ 

OO 

-=t 

LTN 

VO 

r~- 

CO 

ON 

o 

H 

H 
H 

OJ 

rH 

OO 

H 

-d- 

H 

UA 

H 

in 

rH 

OO 

-=f 

O 

OJ 

J- 

CO 

o 

H 

t— 

rH 

OO 

rH 

LTN 

ON 

H 

rH 

vO 

H 

OJ 

H 

•A 

*\ 

^\ 

■N 

*s 

** 

•A 

•s 

»s 

•N 

•\ 

•A 

•N. 

•A 

'A 

*\ 

o 

H 

OJ 

OO 

J- 

LTN 

VO 

t- 

CO 

OA 

O 

rH 

H 

H 

OJ 
H 

OO 

H 

rH 

UA 

H 

LTN 

rH 

OO 

J" 

O 

OJ 

CO 

O 

rH 

c— 

rH 

OO 

H 

LTN 

ON 

H 

H 

VO 

H 

OJ 

rH 

^t 

*\ 

•\ 

»\ 

•N 

•s 

•s 

•a 

»v 

•N 

»s 

»\ 

•\ 

•A 

•A 

•N 

•A 

o 

rH 

OJ 

OO 

-=f 

UA 

VO 

t- 

CO 

ON 

O 
H 

H 
H 

OJ 
rH 

OO 
H 

rH 

UA 

H 

UA 

rH 

OO 

-=r 

O 

OJ 

o 

H 

t- 

rH 

OO 

rH 

UA 

ON 

rH 

rH 

VO 

H 

OJ 

rH 

-Cf 

CO 

»\ 

*\ 

*s 

•s 

»\ 

»\ 

•N 

■N 

•N 

•v 

*\ 

•\ 

•A 

*\ 

•N 

•\ 

o 

rH 

OJ 

OO 

-=t- 

LTN 

VO 

c— 

CO 

ON 

O 

rH 

rH 
rH 

OJ 
H 

OO 
rH 

■H 

UA 
H 

6.   PROGRAM  DESCRIPTIONS 

6.1  Reordering  matrix  elements  into  modified  skew  storage  scheme. 

It  is  assumed  in  this  program  that  the  original  matrix  is  stored  in 
a  straight  storage  scheme  at  the  entry  to  this  program.   The  procedure  named 
DTST  is  therefore  called  for  reordering  elements  in  each  row  of  this  matrix 
(this  procedure  is  a  part  of  this  program).   To  arrange  the  resultant  Hessen- 
"berg  matrix  into  the  straight  storage  scheme,  DTST  is  again  applied.   The  pro- 
cedure DTST  maps  the  jth  element  in  the  row  vector  into  the  jth  position. 
The  algorithm  of  this  mapping  consists  of  the  following  three  steps: 

I.  Exchange  between  elements  whose  binary  expressions  are  0  jYj^jpjn  1  and 
elements  1  j^j-Jo^i  °  for  a11  ccm'bina'fcions  °f  $   ,    .jpJ,  j  ,    and  i  .   This 

is  done  by  the  proper  mode  bit  setting  for  picking  up  the  elements  in  the 
form  of  0  ji  j ojpj-i  -L  arL(^  a  routing  with  distance  31  followed  by  another  mode 
bit  setting  for  picking  up  elements  1  JLJ-JpO-,  0  and  a  routing  with  distance 
33- 

II.  Exchange  between  elements  j  Oj  j  lj   and  Jc-lj^OpOj  .   This  is  done  by  two 
mode  bit  settings  and  two  routings  whose  routing  distances  are  lk   and  50 
respectively. 

III.  Exchange  between  elements  j  j,  Olj  j   and  elements  j  j,  lOj  j  .   This 
step  is  also  accomplished  by  two  mode  bit  settings  and  two  routings  whose 
routing  distances  are  k   and  60  respectively.   The  flow  chart  of  the  procedure 
DTST  is  given  in  Appendix  1. 

After  reordering  elements  of  all  rows  by  applying  DTST,  the  program 
routes  the  i  th  row  of  the  matrix  to  the  right  by  the  distance  i  to  the  matrix 
in  the  modified  skew  storage  form.   It  is  obvious  that  the  unskewing  of  the 
matrix  followed  by  the  application  of  the  procedure  DTST  on  each  row  puts  the 
elements  of  the  matrix  in  the  modified  skew  storage  scheme  into  the  straight 
storage  scheme. 

6.2  Implementation  of  the  Householder's  algorithm. 

N-l     2i 
Mode  bit  setting  for  computing  S  .   The  computation  of  S  =  (  E   a.  J2 

i=r+l  ' 

is   fairly   straight-forward  except   for  the  mode  bit   setting.      The  mode  bit  must 
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be  set  properly  in  order  to  pick  up  a.   for  i  =  r+1, ...,    N-l.  First  we 

1,  r  _ 

look  at  the  content  of  PENBR  at  PE  r  which  contains  r  ,    and  route  PEN  by  the 
distance  r  .   This  row  vector  PEN  is  then  stored  in  register  X(RGX).  We  set 
the  mode  bits  of  all  PE's  whose  register  X  contains  numbers  r+1  through  N-l 


because  these  PE's  contain  a. 


i,r 


at  its  PE  address  A  +  (content  of  register  X) 


The  following  diagram  illustrates  this  step  for  r  =  2,    r  =  16  and  N  =  6k 


PE 

0 

1 

2 

3 

1+ 

PEN: 

0 

1 

2 

3 

h 

PENBR: 

0 

32 

16 

18 

8 

PE  o    1   15   16  17   18   19  20  . . . 

RGX  1+8    h9      63   0   l   2    3  h   ... 

Mode  Bit  1     1   1    0    0    0    1  1  . . . 

PE  Memory 

A:    ao,o  ao,32  ao,6o  ao,2  &o,3k  ao,i8  ao,50  ao,io  " 

al,63  al,0  al,28  ai,6o  al,2  %3h    al,l8  al,50  * 

a2,31  a2,63    a2,¥+      a2,28   a2,60  a2,2        a2,3^  &2,lB' 

a3^T  a3,3l a3,12     *3,Uh  a3,28  a3,6o    a3,2  a3,3^' 

%fl5  %,k-j aU,52    %,12  %,hh  aU,28    %,Go  %,2     ' 


63 

^7 

1 


a 


0,63 
al,31 

a2^7 
a3,15 
%55 


Storage  scheme  used  for  vector  u  and  v 

i  vector  u  ,  w 

The  registers  A  are  cleared  elsewhere.   The  mode 
1  above  and  the  contents  of  PE  address  A  +  (conte 
of  register  X)  are  added  to  the  register  A.   Then  the  register  A  is  routed 


For  finding  the  vector  u  ,   we  place  S  in  the  register  A  of  the  PE 

which  contains  a  , 

r+1,  r 

bits  are  set  as  described  above  and  the  contents  of  PE  address  A  +  (content 
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right  by  6k   -  r  (or  routed  left  by  r)  and  the  contents  of  the  register  A  are 


r,i   2 


stored  in  the  row  named  U.   As  the  result  of  this  calculation,  u^  ..  is  stored 
in  PE  i  in  the  U  vector.   The  vector  V  is  computed  by  multiply] 
the  register  A  and  stored  in  the  row  V  in  the  same  manner  as  U. 


N-l 

Computation  of  P.   The  i  th  element  of  P  ,    P  .is  computed  as  P. =  E    j  ji 

'  j=r+l 

where  i  =  r+1,  .  ..,  N-l.   Since  A  is  stored  in  the  modified  skew  form,  we 

have  to  route  the  vector  U  for  computing  P   . .  We  first  find  i  by  looking  at 

r ,  l 

the  i  th  element  of  PENBR.   The  PE  number  and  vector  U  are  then  brought  into 
the  register  X  and  the  register  S  respectively.   Next  the  contents  of  these 
registers  are  routed  by  i  and  the  mode  bits  of  PE's  which  contain  r  +  1  through 
N  -  1  in  the  register  X  are  set  to  1.   Then  we  multiply  the  content  of  the 
register  S  to  the  content  of  PE  address  (A  +  content  of  RGX)  for  obtaining 
u.a..  for  j  =  r+1,  ....  N-l.   The  sum.  of  all  numbers  in  the  enabled  register  A 
is  formed  and  stored  in  the  i  th  position  of  the  vector  P. 

N-l 

Computation  of  q.   The  i  th  element  of  q,  q   .  is  computed  as  q.  =   E  a.  .  u. 

T*   T";1  i    j=r+1  13   J 

where  i  =  0,  1,  . . . ,  N-l.   Since  A  is  stored  in  the  modified  skew  storage  and 
u   .  is  stored  in  the  straight  storage,  the  following  procedure  is  taken  to 
carry  out  the  computation.   Applying  the  procedure  DTST,  all  elements  of  the 
vector  U  are  reordered  first.   The  i  th  element  of  the  vector  U  is  placed  at 
PE  i  of  the  register  S.   The  vector  PENBR  and  the  register  S  are  routed  by  i  and 
the  former  is  placed  in  the  register  X.   The  mode  bit  of  PE's  whose  register  X 
contains  r+1  through  N-l  are  set  to  1  and  the  product  of  the  register  S  and 
the  content  of  PE  memory  A+i  is  formed  on  the  register  A.   The  sum  of  the 
contents  of  all  enabled  register  A  is  stored  as  the  i  th  element  of  the 
vector  Q,. 
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Appendix  1.   Flow  Chart  of  DTST 
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Appendix  2.  Flow  chart  of  HSBG 
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